Effect of hydrostatic pressure on transport in the topological insulator Bi 2 Te 2 Se 
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The Hall coefficient R H and resistivity p of the topological insulator I^TeiSe display a number of puzzling 
features vs. temperature T . We propose a model that describes well the non-monotonic variation of R H (T). 
A key feature of the model is a chemical potential that is weakly T-dependent. From the fit to the model, we 
infer a "transport" gap Aj- of 50 mV. We find that hydrostatic pressure P (0-27 kbar) has a pronounced effect on 
both R H and p. We show that these changes arise from decreases in both A T and the hole effective mass under 
pressure. 



PACS numbers: 74.62.Fj, 73.20.At, 72.20.-i, 73.25.+i 

Topological Insulators are a class of semiconductors in 
which the bulk energy gap is traversed by current-carrying 
surface states, which have a massless, Dirac dispersion. A 
novel feature of the surface state is the locking of the spin 
of the surface electron perpendicular to its momentum with 
a helicitythat has opposite signs on opposing faces of a 
crystal OJ-la]- The spin locking strongly suppresses back- 
scattering. In several Bi-based semiconductors, the existence 
of the topological surface state (SS) and its spin-locked na- 
ture have been established by angle-resolved photoemission 
spectroscopy (ARPES) i7 Hl0ll . Detailed scanning tunneling 
microscopy (STM) experiments have also confirmed the spin- 
locked nature of the surface states Hill . In transport exper- 
iments, Shubnikov-de Haas (SdH) oscillations from the SS 
have been detected in Bi 2 Te 3 H, and in (Bi,Sb) 2 Se 3 fill . 
Analysis of the SdH oscillations in Bi2Te3 Il4tl yields a sur- 
face mobility p. s = 8,000-10,000 cm 2 /Vs. Despite the high 
mobility, the bulk conductance Gb can exceed the surface con- 
ductance G s by factors of 3000 or more in crystals of these 
materials. A major experimental task is to bring the ratio 
Gb/G s significantly below 1. 

Recently, the hybrid semiconductor Bi 2 Te 2 Se has emerged 
as a highly promising topological insulator (TI) for transport 
experiments lfl6l[l7ll . ARPES measurements fl8ll reveal a sin- 
gle spin-locked Dirac state crossing a bulk energy gap of ~350 
mV. The Dirac point lies close to the top of the bulk valence 
band at the F point. In carefully annealed crystals, the ob- 
served resistivity p at 4 K (Fig. CD) is typically 1 Qcm, but 
can be much larger. In Refs. 1171 [1911 . we show that in crys- 
tals with p(4 K) = 6 Qcm, the surface SdH oscillations are 
especially well-resolved. The surface SdH oscillations yield a 
surface mobility fx s = 2,800-3,200 cm 2 /Vs, which implies that 
the ratio Gb/G s is ~1. At 4 K, the residual bulk carriers are 
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n-type with a bulk density «/, equal to ~ 3 x 10 16 cm 3 (com 
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pared with 10 ls -10 ly cm" 3 in Bi 2 Te 3 and Bi 2 Se 3 ). Hall mea- 
surements show that the bulk mobility pi, is 60 times smaller 
thanp, at4K lfH[l7ll . 

Despite the large p in Bi 2 Te 2 Se, much of its bulk conduc- 
tion features are puzzling. The Hall coefficient Rh displays an 
interesting non-monotonic profile vs. T that seems unrelated 
to thermal activation at first glance (Fig. [2]). To further im- 
prove the quality of Bi 2 Te 2 Se crystals, a better understanding 



FIG. 1: (Color online) The resistivity profiles, p vs. T, at selected 
pressures in Bi2Te2Se for samples SI - • ■ S4 (Panels a- • • d, respec- 
tively). Resistivities of Samples S3 and S4 at higher P are shown in 
Panels (e) and (f), respectively. In sample SI, the curves were mea- 
sured in the sequence P = 0, 1 1.3 and 5.1 kbar. In all panels, p is 
plotted along the vertical axes, and T along the horizontal axes. 



of where the bulk carriers come from is clearly important. We 
propose a model that accounts well for the T dependence of 
the weak-field conductivity tensor. We show that the trans- 
port quantities are consistent with thermal activation across a 
small transport gap of ~50 mV. A key feature of the model is 
a weakly T-dependent chemical potential p(T). Hydrostatic 
pressure allows us to tune the system parameters. We show 
that the observed changes are consistent with the model. In 
addition, we uncover an interesting trend under pressure in 
the magnetoresistance. 

Four crystals of Bi 2 Te 2 Se, grown by a modified Bridgman 
method lflq,[l7ll (Samples SI - ■ • S4), were investigated. Two 
piston-cylinder pressure cells were used, with high-purity Pb 
or Sn serving as the manometer. The pressure was determined 
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FIG. 2: (Color online) The temperature dependence of the Hall coef- 
ficient R H in samples S2 (Panel a) and S3 (Panel b) at the pressures 
indicated. The solid curves in (a) and (b) are guides to the eye. Pan- 
els (c) and (d) plot the conductivity cr xx and Hall conductivity cr„ 
derived from R H and p in S2 at P = (open symbols). The bold 
curves are the fits of the data to the proposed model. In Panel (d), 
we replot for clarity the low-T data (triangles) and the fit (thin curve) 
in expanded scale (xlO). The effective mass m h and transport gap 
A r (0) = Ao derived from the fits in S2 are plotted vs. P in the insets 
in Panels (c) and (d), respectively. 



by the pressure-dependent superconducting transition temper- 
ature of Pb (or Sn), as well as the resistance ratio R{p)/R(0) 
at room temperature. Thin crystals with mirror-like surfaces 
(nominally 1.5 x 0.6 x 0.02 mm in size) were cleaved from 
the ingot. The cleaved faces expose the outer Te layers. The 
axes a and b lie in the layer, and the axis c is normal to the 
cleavage plane. Ohmic contacts were made in the Hall-bar 
configuration on each crystal. With daphne oil as the pressure 
fluid, we attained a maximum hydrostatic pressure of 27 kbar. 

Figure [T]plots the temperature dependence of p under pres- 
sure. At ambient pressure (P = 0), the observed resistiv- 
ity p initially increases steeply as the temperature T is de- 
creased from 300 K, but tends towards a plateau value below 
50 K. This saturation implies a parallel metallic conduction 
channel that has been identified to be the topological surface 
states lfl6l[l7ll . With applied pressure, the p-T profile gener- 
ally becomes much more metallic at all T. We observe two 
classes of behavior. In S2, p at 4 K decreases monotonically 
with increasing P (Fig. Ufa), whereas, in SI, S3 and S4, the 
change is non-monotonic if P < 10 kbar (Panels a, c and d). 
However, above 10 kbar, p(T) tends towards a metallic profile 
with a weak T dependence in all samples. 

Pressure also has a pronounced effect on the Hall coefficient 
R H , as shown in Fig. [2]a (in S2) and Fig. Wp (S3). At P = 0, 
the curve Rh vs. T in S2 displays a characteristic peak at 
80 K and a zero-crossing at 58 K fit)}]. In applied pressure, 
both features shift monotonically to lower T. Simultaneously, 
the large negative plateau at low-T moves upwards, eventually 




FIG. 3: (Color online) Sketch of the band structure near the T point. 
The left panel shows the Dirac surface state (SS) as bold lines. The 
chemical potential ji lies ~130 mV above the Dirac point (DP). 
Shaded regions represent bulk bands. The maximum of the valence 
band lies below /i by the energy spacing A T (the maximum is likely 
at a finite k). The fits show that fj is T dependent. The right panel 
shows the densities of states of the bulk conduction and hole bands 
(N e and N/,, respectively) and the surface states N s . The solid curve 
labelled as f(E - fi) is the Fermi-Dirac distribution. 



attaining a positive value at 11.7 kbar. For S3, Rh is already 
positive at all T at 10 kbar. To exploit the additivity of the 
conductivities, we invert the resistivity tensor p,j to obtain the 
conductivity tensor cry. The inferred conductivity <r xx and 
Hall conductivity cr^ are plotted as open circles in Panels (c) 
and (d), respectively. Also plotted are the fits (bold curves) to 
these quantities using the following model. 

As may be seen in the expanded scale in Panel (d), cr xy 
(open triangles) is negative and nearly T-independent below 
30 K. Above 30 K, cr xv undergoes an exponential increase to 
very large positive values. This suggests that we have a small 
density of n-type carriers at low T, which are overwhelmed by 
a divergent population of holes induced by thermal activation 
across a small gap. However, because cr Ay rises to a broad 
peak at 180 K, the exponential growth is eventually tempered 
by high-T processes. 



TABLE I: Fit parameters at selected pressures P (for Sample S2). Ao 
is the transport gap at T = 0, y the coefficient of the T 2 term in the 
gap, mi, the hole effective mass (with mo the free electron mass) and 
P is the exponent in the hole mobility. 
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At 4 K, we assume that charge conduction involves carriers 
in the topological surface state and a dilute concentration of 
bulk carriers (both n-type). As found in Ref. fl7il . the bulk 
carriers, which have very low mobility (pi, ~ 50 cm 2 /Vs), oc- 
cupy an impurity band. By contrast, the high mobility of the 
surface electrons (p s ~ 3,000 cm 2 /Vs) leads to SdH oscilla- 
tions (see below). The data at 4 K, restricted to relatively low 
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B, do not allow the surface and bulk electrons to be separated 
(in Refs. lfl4l[l7tl . the separation was accomplished by using 
the SdH oscillation amplitudes), we treat the two bands here 
as one band with an effective mobility fi e = A/(l +CT a ), with 
a, A and C as free parameters. 
The conductivity is cr xx 

= 2e 2 2 k (-d/7<9£)v A (k) 2 T, where 
f(E~n) is the Fermi-Dirac distribution, v(k) the group veloc- 
ity at the wave vector k and t is the transport lifetime. With 
the chemical potential /i in the bulk gap as sketched in Fig. 
[3j we obtain for the holes cr xx = nyfifih and <j xy = n^efX^B, 
with /J/, the hole mobility. In the limit A T /k B T » 1 (k B is 
the Boltzmann constant) the hole population has the activated 
form 



1 / 2mi,kBT\ 2 



"" (r) =4\ ntt 



,-& T /k B T 



(1) 



where m/, is the hole effective mass. The energy scale A7- = 
fi - Eb is the "transport gap" that dictates the exponential 
growth in cr xv at low T (Eb is the energy at the valence band 
maximum). As we show below, we may ignore thermal acti- 
vation to the conduction band for T < 300 K. We express the 
hole mobility as /y/, = DT~&, with D and B as free parameters 
(the impurity scattering rate for the holes plays no role in the 
fitting because the holes are thermally activated). 

To fit the data at high T, however, we further assume that 
A T has the temperature dependence 



A T (T) = A + yT 2 



(2) 



with Ao and y as free parameters. 

The fits to the conductivity tensor are plotted as bold curves 
in Figs. [3J: and [2jl. We remark that our primary goal is to 
describe the weak-field hole conductivity tensor using the 4 
constants Ao, y, mi % and B. The strong thermal activation at 
low T and the pronounced non-monotonicity of cr xv at high 
T severely restrict the physically reasonable values that these 
parameters may assume. Hence the convergence of the fits 
is rapid and relatively unambiguous (by contrast, the n-band 
parameters a, A and C are less reliably obtained because the n- 
type carriers are resolved only below 30 K). We have extended 
the fits to finite P as well for S2. The optimal parameters are 
reported in Table I. 

From the fits, we infer that the steep exponential increase 
in cr xv when T is raised above 30 K (triangles in Fig. |2}l) 
results from the thermal activation of holes into the valence 
band (Eq. [TJ. The monotonic decrease of the hole mobility /i/, 
with increasing T tends to counter this increase, especially in 
<T xy ~ [il. However, to produce the non-monotonic variation 
in cr xv , we need A7- to increase slightly with T as in Eq. [2] 
We attribute this increase to a T-dependent chemical poten- 
tial. Particle number conservation dictates that fi must change 
with T when k B T is large enough to reach a band with a large 
density of states (here the hole band). This causes /i to move 
away from the valence band with an increment varying as T 2 . 
At ambient P, the increment is 6.3 mV when T = 100 K (Table 
I). As a consequence, «/,(T) grows much more slowly above 
200 K, resulting in an overall decrease in <x vv as observed. We 
expect a T-dependent p: to be a general feature in topological 
insulators because p lies close to one of the bulk band extrema. 
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FIG. 4: (Color online) Comparison of the magnetoresistance curves 
at ambient pressure and high pressure. In Panel (a), the derivative 
dpxx/dB at P = and T = 0.35 K displays SdH oscillations from 
the surface states. Panel (b) shows that the oscillations are periodic 
in with a period that yields k F = 0.040 A~'. The oscillations 
become unresolved at 11.7 kbar. In Panels (c) and (d), we plot the 
MR curves (measured with B||c) at selected pressures for Samples 
S2 (at 4.6 K) and S4 (2 K), respectively. In both samples, the MR 
changes from a saturation trend (at low P) to one that is large and 
quadratic in B at large P. 



As in Refs. ITil fl6l [l7tl . we may estimate the position of 
by detecting the weak surface SdH oscillations (Fig. |4^). 
From the period of the oscillations at 0.35 K, we obtain the 
surface Fermi wavevector kf = 0.040 A -1 . Using the surface 
Fermi velocity Vf = 6xl0 5 m/s U3.[H], the kp value implies 
that fi lies 130 mV above the Dirac point. ARPES UK] shows 
that the direct gap at the F point (k = 0) is 350 mV. Hence 
the conduction band minimum lies 220 mV above //. As this 
separation exceeds A7- by a factor of ~4, we are justified in 
ignoring thermal excitations into the conduction band. In ad- 
dition, we note that A7- is quite small compared to the spac- 
ing of /i from the Dirac point (130 mV). This implies that the 
valence-band maximum is at a finite k, as drawn in Fig. [3] 

From the fits, we now see that the dramatic effect of P on 
the transport properties stems primarily from the steep de- 
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crease of the transport gap. In Sample S2 (inset in Fig. |2}l), 
the zero-Kelvin value Ao decreases from 5 1 mV (at P = 0) to 
15 mV (P = 1 1 .7 kbar). Thus, at 1 1 .7 kbar, the smallness of A 
results in a large population of excited holes even at 4 K (Rh 
is close to zero at 1 1.7 kbar). Nonetheless, a vestige of the ac- 
tivated behavior can be seen at 1 1.7 kbar in p vs. T (Fig. [TJ:). 
At larger P (16-27 kbar), At is nearly completely suppressed 
(Fig. QJl, for S4), and the system remains p type down to 4 K. 
In addition, the fits reveal that an increasing P decreases the 
hole effective mass from 0.16 times the free mass (at P — 0) 
to 0. 1 1 (at 11.7 kbar). The most significant effect, however, is 
that increasing P suppresses the surface conductance, render- 
ing it difficult to resolve against the increased hole conduction. 
Direct evidence for this suppression is obtained by comparing 
the curves of dp xx /dB (Fig. Ek,b). The prominent SdH os- 
cillations at P — become unresolved at 11.3 kbar. Thus, in 
closing the transport gap, pressure converts the system from a 
TI in which surface conductance G s is easily observed at 0.35 
K to one that is dominated by the bulk holes. 

Tuning the relative weights of G s and Gb by pressure pro- 
vides a powerful way to separate transport characteristics of 
the surface states from bulk states. As an example, we exam- 
ine the transverse magnetoresistance, defined as the fractional 
change in p(B), viz. MR{B) = [p(B)/p(0)] - 1, measured 
with B||c. Figures @t andHJl plot the curves of MR(B) at sev- 
eral pressures for Samples S2 (at 4.6 K) and S4 (at 2 K), re- 
spectively. In both samples, when G s and Gb are comparable 
(at low P), the MR tends towards saturation at large fields. 
However, when G/, dominates (large P), the MR switches to 
a quadratic variation (~ B 2 ) consistent with the semiclassical 



MR observed for bulk states in a high-mobility semi-metal. In 
both samples, the dominance of the bulk conduction at large 
P is again evident in the large, semi-classical B 2 variation of 
p. Currently, there is very little understanding of the MR of 
both the bulk and surface states, as discussed in Ref. fTill . 
The present results suggest that the behavior of the MR with 
pressure tuning may reveal specific characteristics of the topo- 
logical surface state, and provide a way to readily distinguish 
surface from bulk conduction. 

Finally, we remark that the fits provide a clearer picture of 
the overall band structure in crystals of Bi2Te2Se that have 
p inside the bulk gap. We obtain an accurate determination 
of the gap Ar(T), which dictates the bulk carrier concentra- 
tion above 30 K, as well as its sensitivity to applied pressure. 
It should be instructive to see how these numbers are altered 
when the pure compound is doped with donors or acceptors. 
An interesting test of the inferred band alignments may be 
made with thermopower, which is sensitive to the band pa- 
rameters. 
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